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SENSITIVITY ANALYSIS OF TORSIONAL VIBRATION 
CHARACTERISTICS OF HELICOPTER ROTOR BLADES 

Part II. AERODYNAMICS AND SENSITIVITY ANALYSIS 

By Theodore Bratanow* and Akin Ecer** 
University of Wisconsin-Milwaukee 


SUMMARY 


A theoretical investigation of dynamic response characteristics of ro- 
tor blades was carried out with special emphasis on torsional degrees-of- 
freedom. The major steps of the analysis are; 

determining the aerodynamic loads on the blade for selected 
forward flight conditions 

. application of a numerical method for determining the sen- 
sitivity of blade vibration characteristics with respect to 
structural, geometric and aerodynamic properties and flight 
conditions. 

Coupled equations of motion for flapwise bending and torsion were for- 
mulated at varying azimuth positions for rotor blades with noncollinear aero- 
dynamic, clastic and mass axes. Both structural and aerodynamic mass, damp- 
ing and stiffness coefficients were included. The vibrations of a’ sample 
blade at different flight conditions were investigated from these equations. 
The obtained numerical results were illustrated. 

The sensitivity of overall blade vibration characteristics to torsional 
oscillations was also investigated from the equations of motion for the 
sample blade. The illustrated results show the importance of torsional de- 
grees-of-freedom in rotor blade analysis. Various possibilities of improv- 
ing the overall response by tuning blade geometric, structural and aerody- 
namic characteristics are discussed. 


* Professor, ** Assistant Professor, Department of Mechanics, College of 

Engineering and Applied Science. 


INTRODUCTION 


The objectives of the conducted investigation were to determine the var- 
iation of torsional oscillations and their effect in forward flight on blade 
dynamic response characteristics and to analyze the sensitivity of torsional 
oscillations to blade aerodynamic, structural and geometric properties. In 
the first report [1] , the structural aspects of the analysis were presented. 
In this report the aerodynamic characteristics of rotor blades are defined 
for varying forward flight conditions. The unsteady effects due to the 
pitching motion of the blade were determined from two-dimensional airfoil 
theory [2] . The aerodynamic loading in each of the structural vibration 
modes was determined for different flight conditions and then the aerodynamic 
stiffness and damping coefficients were calculated. These coefficients were 
compared with structural stiffness and damping coefficients for varying azi- 
muth positions. Then the structural and aerodynamic coefficients were com- 
bined for a complete, coupled vibration analysis in forward flight. The sen- 
sitivity of the torsional vibrations in forward flight to structural and geo- 
metric properties of the sample blade was investigated. The effects of the 
variations of mass, elastic and aerodynamic axes of the sample blade were il- 
lustrated for various flight conditions. 

SYMBOLS 

a lift curve slope for forward flow 

a* lift curve slope for reverse flow ’ 

a^j generalized aerodynamic mass coefficient 

b^j -generalized aerodynamic damping coefficient 

B tip loss factor 

blade chord-width, m. 

generalized aerodynamic stiffness coefficient 
lift coefficient 

generalized control moment stiffness in torsion 
damping matrix 

f^Ct) generalized external load function in each mode, N 
generalized external load function, N 
generalized force in bending, N 
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generalized control moment in torsion, N^ih ■ ' ' 

h vertical displacement of the blade at the elastic center, mi 

h. eigenfunction corresponding to- the bending displacement component 
of 2^, m. 

stiffness matrix 

L lift generated by a blade, N 

lift generated by a blade in reverse flow, N 
L . forcing function component for lift, N 

" 2 

m^j^ generalized coupled torsional stiffness in bending, N-m 
Mei(iJj) forcing function component for moment, N-m 
Mp a feathering moment, N-m 

M mass matrix ' ' - ' ' ' 

eigenvector of the- structural system ' ' 



P. . 
ij 



r^ 

R 


R 


ki 



U 

P 


U 


t 


generalized torsional damping' in tors ion, N-sec./rad • 

generalized bending damping in bending, N-sec./m - 

generalized torsional stiffness in torsion, • N/rad 

generalized bending stiffness in bending, N/m 

distance of rotor blade element from center of rotation, m, 

coupled displacement vector, m. 

rotor radius, ra. 

•I 

generalized coupled bending stiffness in torsion, N/m 
generalized coupled bending damping in torsion, N-sec./rad 
time Csec.) 

normal component of the resultant velocity at a blade element, m/sec. 

tangential component of the resultant velocity at a blade element, 
m/sec. 
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induced velocity, m/sec. 

V forward velocity, m/sec. 

Xq eccentricity between the elastic and mass center, m. 

u 

generalized bending deflection at blade tip, m. 

generalized torsional rotation at blade tip, rad 

generalized displacement in a coupled mode, m. 

generalized structural mass coefficient 

blade element angle of attack, rad 

3 blade flapping angle at particular azimuth position, rad 

3- eigenfunction corresponding to the bending rotation component of , 

^ rad ^ 

3^ generalized structural damping coefficient, N-sec./m 

generalized structural stiffness coefficient, N/m 

e an angle defining reverse flow region, rad 

0 pitch change due to blade torsion, rad 

0^ lateral cyclic pitch control input, rad 

0. eigenfunctions corresponding to the torsional rotation component of 

^ Ei’ 

0^ collective pitch, rad 

0^ longitudinal cyclic pitch control input, rad 

0p pitch angle of a blade element without considering torsion, rad 
0j built-in blade twist, rad 

Kc distance between elastic axis and center of pressure in reversed flow 

region, m. 

X inflow ratio (see Eq. (A6)) 

y tip-speed ratio 

q ratio of torsional rotations to bending displacements in the vibration 

mode corresponding to the lowest frequency, rad/m 
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3 

p mass density of air, kg/m 

azimuth angle, rad 

0) coupled frequency of the blade, rad/sec. 

n angular velocity of a rotating blade, rad/sec. 

Subscripts 
A aerodynamic 

s structural 

Superscripts 
B bending 

t transpose 


T 


torsion 



AERODYNAMICS OF ROTOR BLADES IN FORWARD FLIGHT 

Background 

The effects of torsional oscillations of blades become an important con- 
sideration in predicting dynamic response characteristics of high performance 
rotorcraft. Torsional oscillations of rotor blades have been considered by 
researchers in connection with stability and stall problems. The stability 
of coupled bending and torsional vibratory motion of rotor blades at high 
forward speeds has been treated by several investigators [3,4,5]. The tor- 
sional loads used in such analyses have been calculated at lower tip-speed 
ratios, when the coupled torsional rotations are relatively small. Under 
these conditions and when the blade has a sufficiently high torsional stiff- 
ness compared to the bending stiffness, the torsional motion has a minor ef- 
fect on blade stability. Large torsional rotations change the aerodynamic 
characteristics considerably and affect the stability of the rotor blade. 

The stall characteristics of rotor blades have also been studied [6,7, 
8]. Ham [7,8] for instance, has shown that for stall conditions the aero- 
dynamic damping of torsional oscillations becomes negative as the pitch 
angle increases. To balance such a negative aerodynamic damping, he sug- 
gested measures for addition of a positive mechanical damping. 

Various investigators have been involved in determining bending and 
torsional aerodynamic loads on vibrating blades in forward flight [9,10]. A 
complete investigation of the response of rotor blades at high forward speeds 
requires a three-dimensional aeroelastic representation of the problem. The 
aerodynamics of a three-dimensional flow around the rotor blade at high for- 
ward speeds is not sufficiently developed for an accurate determination of 
coupled aerodynamic loads. 

In the conducted investigation, effects of torsional oscillations were 
evaluated by quantitative comparisons of a series of structural and aero- 
dynamic parameters. The aerodynamic loading in forward flight was deter- 
mined according to the classical rotor blade analysis [3,4,11,12,13] as de- 
scribed in the Appendix. 


Method of Solution 

The equations of motion of a rotor blade can be expressed in a general 
matrix form as follows: 


Mr + Dr + Kr = F (t) 


( 1 ) 


where the mass, damping, and stiffness matrices can be considered as con- 
sisting of two components, structural and aerodynamic. 
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M = + 4(4') (2) 

K = 4 + 4(4') (3) 

D = 4 + 4(4') (4) 

The aerodynamic terms in Eqs. (2), (3) and (4) are functions of blade azi- 
muth position and flight conditions. Considering only the structural com- 
ponents of the matrices in Eqs. (2), (3), and (4) and generating a struc- 
tural damping matrix proportional to the critical damping, one can obtain 
the eigenvectors for the following system: 

2 

-03 Mp. +Kp. =0 where r = Z.p. (5) 

— S'^l — S*-l — 


so that the generalized structural mass, stiffness, and damping terms can 
be calculated as 


and 


a. = p. M p. 
1 — s 

3. = £.- 4 Ei 


( 6 ) 

(7) 

( 8 ) 


p^^M p. =p^D p.=p^K p. = 0 for i X j 


(9) 


The generalized aerodynamic mass, damping, and stiffness terms are defined 
as 


= E* tU % 
bij W = eI 5a Ej 

Cij W ' Ei 5a Ej 


(10) 

( 11 ) 

( 12 ), 


Then one can write the generalized equations of motion as 


“i " ^ ^i M ^i " } h* ] 4 " } h = 


11 


(13) 
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In Eq. (13), the aerodynamic terms are coupled in each mode. The response 
problem can be analyzed from the solution of the system of second order dif- 
ferential equations with periodic coefficients in Eq. (13). These coeffi- 
cients consist of a constant uncoupled term and a periodic coupled term. 

The aerodynamic coefficients in Eq. (13) were determined for the flow 
regions shown in figure 1 in terms of the coupled vibration modes as de- 
scribed in the Appendix. These coefficients were calculated for typical 
flight conditions. 


Analysis of the Aerodynamic Loading in Bending and Torsion 

A sample blade (CH-34) was chosen to demonstrate the variation of the 
aerodynamic coefficients for different modes. The structural generalized 
mass and stiffness terms for each mode and coupling between modes were cal- 
culated as shown in table 1. The structural coupling between the lowest 
three modes is summarized in table 2. To illustrate the aerodynamic cou- 
pling between each of these modes, the aerodynamic stiffness and damping 
coefficients for the different modes were calculated as shown in figures 2, 
3, and 4. From the curves one may observe the degree of coupling between 
the aerodynamic forces in each mode. 

Each of the terms defined in Eqs. (A19) , (A20), and (A21), was plotted 
for various tip-speed ratios for the mode corresponding to the lowest bend- 
ing frequency. The aerodynamic bending stiffness, represents the re- 
sultant force in the first structural mode due to unit displacement in the 

same mode and is plotted in figure 5. For a tip-speed ratio corresponding 

to the indicated operational speed of a CH-34 blade (y = 0.22), the aero- 
dynamic bending stiffness of the blade is in the range of the calculated 
structural stiffness. For this case the effect of the reverse flow region 
is negligible. In general, as the tip-speed ratio increases, the bending 
stiffness values show a greater variation with changes in azimuth positions; 
this variation of the stiffness is greater in the forward flow region and 
smaller in the reverse flow region. Therefore, for high tip-speed ratios of 
the blade, stability problems may arise in the flapping mode due to the re- 
sulting negative bending stiffness. 

As shown in figure 6, the generalized aerodynamic damping in the first 

bending mode, also depends on the tip-speed ratio of the blade. For 

the operating speed of the CH-34 blade, the aerodynamic damping coefficient 
is found to have nearly a constant value; however, it varies consid- 
erably when the tip-speed ratio is increased. The aerodynamic damping in- 
creases in the forward flow region but it nearly vanishes for the reverse 
flow region. The decrease in the aerodynamic damping may be one of the 
reasons for the existence of increased bending amplitudes in the reverse 
flow region. It should be noted also that the calculated values of aero- 
dynamic damping are found to be considerably high compared to the struc- 
tural damping of the CH-34 blade. Such results suggest that a further 
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examination of values of characteristic aerodynamic damping for rotor blades 
at high forward speeds is desirable. 

The aerodynamic stiffness term from the coupling of the first bending 
and torsional modes was also found to be dependent on the forward speed. As 
shown in figure 1 , this stiffness term has a different sign for different 
flow regions at high tip-speed ratios, but it is nearly constant for low tip 
speed ratios. 

The variation of torsional stiffness and damping parameters with azi- 
muth positions and tip-speed ratios is shown in figures 8, 9, 10, and 11. 
Aerodynamic torsional stiffness in the first mode exists only in the reverse 
flow region. As can be seen from figure 8, this stiffness acts as a neg- 
ative feathering spring for the reverse flow region and its magnitude is pro 
portional to the increase in tip-speed ratios. Accordingly, the total tor- 
sional stiffness of the blade has lower values in the reverse flow region. 
The structural torsional stiffness of the blade must be increased in order 
to overcome this negative aerodynamic stiffness. 

The variation of the torsional damping for various tip-speed ratios 
was also investigated and is shown in figure 9. For high tip-speed ratios, 
torsional damping becomes higher for the forward flow region, while for the 
reverse flow region the aerodynamic torsional damping vanishes. The vari- 
ations of the coupled aerodynamic stiffness and damping terms are shown in 
figures 10 and 11. These coupled terms also exist only for the reverse 
flow region and increase proportionally to the tip-speed ratios. 


ANALYSIS OF THE COUPLED EQUATIONS OF MOTION 
OF A ROTOR BLADE IN FORWARD FLIGHT 

Method of Solution 


The coupled equations of motion of rotor blades in bending and torsion 
can be expressed in terms of the generalized mass, damping and stiffness 
coefficients. These equations are: 


X. 

1 1 


+ 

3? + P. -HlJ) 

X. . + 

Y? + Q. . (il^) 1 


H 

H 

__1 

1.1 

H 

H 


X. . + Y, m.,.(ij;) = L . (ijj) 

( 14 ) 


for the bending motion, and 



T •• 

a. Y. + 

bJ + p.,(>i^) 

Y. . + 

yJ + q. . (i|^) 

1 1 

1 ij 

1,3 

1 1 j 






+ SihW Xh = MeiW ■ (15) 

for the torsional motion where (X^, Y^) is the coupled structural mode. 

Eqs. (14) and (15) can be treated as a system of differential equa- , 
tions with periodic coefficients to analyze the dynamic response of rotor 
blades at high forward speeds. An approximation of the response of the 
coupled system can be examined from the analysis of a two-degree-of- 
freedom system. 

The displacements and rotations of the blade can now be assumed as 
follows : 

X = X e“^ (16) 

1 o ^ 

and 

Y = Y e“''^ = nX (17) 

1 o ' o 

where n is the ratio of torsional rotations to bending displacements. Eqs 
(14) and (15) can be written for the free vibrations as follows: 


2 B 

0) + 




p^tw 


0) + 


^1 " 


QnW 


nm^^(t|;) = 0 


(18) 



T 
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(jJ + 


T 

^1 " 


qii(^) 


+ Rj^(ip) + S^^(iJ;) to = 0 

(19) 


These two algebraic equations can be solved to obtain the following ex- 
pressions: • 


n 


2 T 

to + 


R^^(t|j) + to Sjj(i|;) 

^1 Pll^'*^^ w + Y j + q^j(\jj) 


( 20 ) 


and 


10 


0) 



Pllw 


B 

w + Y, + 


Qll(« 


2 T 

iti + 






'll 







( 21 ) 


Eq. (21) is a fourth order polynomial which can be solved to obtain the 
corresponding lowest frequencies of the coupled equations of motion. From 
Eq. (20) the degree of coupling can be calculated. One has to consider the 
variation of each one of the coefficients of Eqs. (20) and (21) in order to 
explain the obtained results. According to the azimuth position of the 
blade, the coefficients of Eqs. (20) and (21) exhibit different character- 
istics. 

In forward flow region - For zero aerodynamic damping (k = 0.5) in the 
forward flow region, the coupled damping and stiffness terms and S^^^) 

vanish as shown in figures 9 and 10. Then, Eqs. (18) and (19) become uncou- 
pled, as can be seen from Eq. (21). The torsional oscillations can be de- 
termined for this case only from the structural and aerodynamic properties 
of the blade in torsion. However, for an eccentricity between elastic and 
aerodynamic axes, (k / 0.5), a complete solution of Eq. (21) is necessary. 

In reverse flow region - When the blade enters the reverse flow region, 
the magnitudes of the coupled terms are considerably higher. The coupled 
motion of the blade can then be analyzed from the complete solution of Eq. 
(21) . The coupled frequencies are strongly dependent on the degree of cou- 
pling in this region. 

The magnitudes of the torsional oscillations can be analyzed in terms 
of the following two factors; 

2 2 

a) the ratio of torsional to bending frequency (cOr^/Wg) 

b) the degree of coupling between torsional and bending modes. 

If the magnitudes of the torsional and bending frequencies are close, 
a major portion of the torsional oscillations occur in the mode corre- 
sponding to the. torsional frequency. As can be seen from figures 7 and 8, 
the value of the torsional frequency 0 ).p depends on the negative aerodynamic 

stiffness and the decrease in the aerodynamic damping for the reverse flow 
region. When the torsional frequency is much higher than the bending fre- 
quency and the coupling between bending and torsional modes is significant, 
the torsional oscillations are mostly in the coupled mode corresponding to 
the lowest bending frequency. The presence of torsional vibrations, cou- 
pled with these bending oscillations, is indicated from the results of the 
harmonic analysis and from experimental measurements of pitching and 
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pliinging motion of rotor blades [8]. 


Sensitivity Analysis of Torsional Vibration Characteristics 

of a Sample Blade 

The sensitivity analysis of coupled torsional and bending vibrations was 
carried out for a sample blade (CH-34) as described above. Quantitative re- 
sults were obtained from the solution of Eqs. (20) and (21) for different 
flight conditions and varying the blade properties. 

Figures 12-14 show the response of the blade in forward flight for coin- 
ciding elastic and aerodynamic axes. For the operational speed of the 
sample blade (U = 0.22), the aerodynamic damping and the natural bending 
frequency show a periodic variation, corresponding to the variation of bend- 
ing stiffness and damping coefficients in figures 2 and 3. As shown in fig- 
ure 9, torsional damping varies also periodically with changing azimuth po- 
sitions. However, the constant torsional structural stiffness is the im- 
portant factor in determining torsional frequency, as can be seen from fig- 
ure 8. The coupling ratio between torsional and bending modes (q) exists 
only in the reverse flow region and reaches a maximum at an azimuth position 
of 270°. 

A detailed torsional sensitivity analysis of the rotor blade was further 
developed starting from the ideal uncoupled case. The first term to be in- 
cluded in the sensitivity analysis was the effect of the eccentricity be- 
tween elastic and mass axes. The effects of structural coupling bn blade 
response, considering only the structural terms, were discussed in the first 
report [1]. The analysis can be extended to include the additional aerody- 
namic uncoupled loading of the blade. As shown in figures 15-17, in this 
case, bending and torsional frequencies were insensitive to the addition of 
structural eccentricity of the blade. Such behavior can be explained in 
terms of the relative magnitudes of structural and aerodynamic coefficients. 
As described in the first report, bending frequencies are not sensitive to 
the structural eccentricity (xg) even without considering the aerodynamic 

loading. The importance of this effect is indicated in the variation of 
coupling ratio with varying eccentricities. The coupling in the bending 
mode is increased uniformly for all positions of the blade for an eccentric- 
ity of Xg = -1, and it exists mostly in the forward region for Xg = 1. The 

coupling ratio for the torsional mode is increased with additional eccen- 
tricity. The torsional frequencies are also increased as shown in figure 15 

The second term considered in the torsional sensitivity analysis of the 
blade was the effect of the relative position of aerodynamic axis with re- 
spect to the elastic axis. The exact position of the aerodynamic axis var- 
ies with time as the azimuth position of the blade changes in forward flight 
The results obtained from the sensitivity analysis can be used in this case 
to evaluate the effects of the changes in the aerodynamic and elastic axes, 
both qualitatively and quantitatively. The previous results were obtained 
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for an ideal case, when the blade oscillates around the one-quarter chord- 
width and the distance between the elastic and aerodynamic axes is one-half 
chord-width for both reverse and forward flows (k = 0.5). To investigate 
the effect of the changes in the relative positions of the aerodynamic and 
elastic axes, the aerodynamic axis was moved in a symmetric fashion for for- 
ward and reverse flow by changing k to 0.4 and 0.6 respectively. The addi- 
tional aerodynamic pitching moments were included in the expressions in the 
Appendix for varying positions of elastic and aerodynamic axes. The ob- 
tained results were illustrated in figures 18-21. 

Two cases (k = 0.6, k = 0.4) describe the behavior when the aerodynamic 
axis is fore and aft of the elastic axis. The bending frequencies of the 
rotor blade are reduced considerably for the entire rotor disc when K = 0.6. 
However, for K = 0.4 these frequencies show a consistent variation for dif- 
ferent azimuth positions. Compared to the case K = 0.5, the response of 
the blade is stabilized for azimuth positions 90° < < 270° for ic = 0.4. 

This improvement is due to the additional stiffness obtained from the cou- 
pling with the torsional mode. On the other hand, the coupling ratio has 
been increased considerably compared to the uncoupled case (k = 0.5). For 
K = 0.4, the blade also exhibits an underdamped behavior for azimuth posi- 
tions 90° < < 360°. 

The variation of torsional frequencies with coupling between mass and 
elastic axes is also an important consideration in the response analysis. 
While torsional frequencies are increased for k = 0.4, the blade is un- 
stable in a basically torsional mode for.K = 0.6. These results indi- 
cate that the blade vibrations may become unstable in a predominantly tor- 
sional mode, when the aerodynamic axis moves ahead of the elastic axis.. 

The torsional vibrations are overdaraped for k = 0.4. 

From the above results, the sensitivity of torsional characteristics 
of the CH-34 blade was described for the blade operating speed (y = 0.22). 
Another important element of the investigation of torsional sensitivity 
characteristics of rotor blades is the changing flight conditions. The 
torsional and bending vibration frequencies and the coupling between the 
modes were calculated for different flight conditions and different po- 
sitions of the aerodynamic and elastic axes. Figures 22-27 indicate the 
changes in the bending frequencies and coupling ratios with increasing 
forward speeds. 

As the tip-speed ratio increases, the bending frequencies show a peri- 
odic variation with changing azimuth positions and the values of the fre- 
quencies are highest at the forward flow region. These results can be 
interpreted again in terms of the values of bending stiffness in figure 2. 
The bending vibrations become unstable for azimuth position 90° < ijj < 180°. 
The frequencies in the reverse flow region remain relatively low. For 
K = 0.5 the coupling ratio exists again only for the reverse flow region. 
Another possibility of instability for 270° < < 360° exists in the re- 

verse flow region due to the low torsional frequencies in this region. 

For K = 0.4 the effect of tip-speed ratios is considerably smaller in 
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the bending modes due to the aerodynamic coupling terms in bending and tor- 
sion. In this case coupling improves the blade response by removing the in- 
stability region. In the reverse flow region the frequencies show a more 
rapid change as compared to K = 0.5. 

For K = 0.6 the blade becomes unstable in a mainly torsional mode. The 
bending frequencies vary rather strongly as the tip-speed ratios increase, 
due to the coupling with the unstable torsional motion. The ratio of cou- 
pling shows a similar unstable behavior. 

Variations of torsional frequencies and coupling ratios of rotor blades 
for different forward flight conditions are illustrated in figures 28-33. 

In the reverse flow region, the torsional frequencies change considerably 
with increasing tip-speed ratios'. The instability in the reverse flow re- 
gion can be explained in terms of the negative stiffness coefficient in fig- 
ure 8. The coupling ratio for the torsional displacement indicates an un- 
derdamped behavior and varies slightly with increasing tip-speed ratios. 
Torsional frequencies are found to be sensitive to tip-speed ratios, also 
for K = 0.4. The coupling increases the torsional frequencies in the re- 
verse flow region compared to the uncoupled case K = 0.5. The coupling ra- 
tio for this case indicates an underdamped behavior for high tip-speed ra- 
tios . 

For K = 0.6 the effect of tip-speed ratios was also tested. The tor- 
sional frequencies decrease proportionally to the increase in the tip-speed 
ratios. The ratio of coupling for this case indicates an underdamped behav- 
ior and shows a change in the sign of the coupling ratio in forward and re- 
verse flow regions. 

After discussing the effects of structural and aerodynamic coupling 
separately, a combination of eccentricities between elastic, mass and aero- 
dynamic axes is presented as shown in figures 34 and 35. In this case, for 
K = 0.4, y = 0.22, the structural eccentricity Xg was varied for Xg = -1, 0, 

1. As can be seen from figure 34, the bending and torsional frequencies 
were affected only slightly with this structural eccentricity. Such behav- 
ior can be explained in terms of the sensitivity of rotor blades to struc- 
tural eccentricities, as described in the first report, and the relative 
magnitudes of structural and aerodynamic loads, as discussed in this report. 


Conclusions 

The presented investigation of coupled bending and torsional vibration 
characteristics of rotor blades has shown the effectiveness of sensitivity 
analysis for studying dynamic response problems. The importance of the tor- 
sional degrees-of-freedom was established in terms of structural and aerody- 
namic components of the aeroelastic analysis of rotor blades. 

The sample blade analysis has shown that the aerodynamic and structural 
coupling terms in the equations of motion have to be considered to obtain 
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realistic results for vibration frequencies and mode shapes. From the de- 
veloped sensitivity analysis of the mathematical model the significance of 
blade design parameters can be determined and used to improve the blade 
response and stability. The investigation has also shown the importance of 
the aerodynamic terms in the vibration analysis of blades in forward and re- 
verse flow regions at high tip-speed ratios. 


15 



APPENDIX 


CALCULATION OF THE AERODYNAMIC LOADING IN FORWARD FLIGHT 


The airflow over the entire rotor disk, 
analyzed in the three regions shown in figure 
blade encounters only forward flow; in region 
the mixed flow region (2), blade elements for 
to r/R = -y sin ip operate in the reverse flow 
to r/R = 1 operate in normal flow. The angle 
region (2) is defined as 


e = 


/ 

A ■ B 

Arc sin — 

< 

]L 

2 

V 


at any azimuth position, was 
1 [13] . In region (1) the 
(3) only reversed flow. In 
radial positions from r/R = 0 
and those from r/R = -p sin 
defining azimuth positions in 


B < p 


B > p 


(Al) 


Blade Flapping Motion : - From the two-dimensional theory the lift co- 
efficient of a blade element is assigned to be [14] 


where is the angle of attack given 

U 

a = 7^+6 + 

r P 


(A2) 


(A3) 


U 


^ and Up are the tangential and normal components of the resultant air ve- 
locity respectively, and they can be written as 


U^ = f^R + pflR sin ip 


(A4) 


U = AJIR + ^ - pfiR3 cos \p 
p dt ^ 


(A5) 


where 


X = 


V sin a 
r 

ftR 


V. 

1 


(A6) 


and 


_ V_ 
^ flR 


(A7) 
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The pitch angle of the blade section is defined by 


9t,= 0d + 0 = 0 + 0 sin ip + 6 cos ^ 0. +0 (A8) 

Tr os c Rl 

where 0p represents the angle of pitch control and 0 is the pitch angle due 
to torsion. 

The torsional rotation 0, the bending displacement h and rotation 3 
can be defined in terms of the uncoupled structural modes of a rotating blade 
as follows: 


0 = 

Y.0. 
1 1 

(A9) 

h = 

X.h. 
1 1 

(AlO) 

3 = 

x,3, 

(All) 

where 0 . , h . , and 3 ■ are the 

eigenfunctions obtained from Eq. C5). 

For two- 


dimensional flow the lift on the blade element can be written as 

dL = y p c dr (A12) 

where V is the resultant velocity at the blade element. Since it is as- 
sumed that U /U. is small, 

p t 

V = (A13) 


The generalized forces in bending can be written, using the principle of 
virtual ’work, as 


or 





c h dr 


(A14) 



h. dL = 
1 



^ p Cj c h. dr 

^ Ij 1 


(A15) 


Substituting Eqs. (A2) and (A13) into Eq. (A15), the generalized lift force 
in each mode can be written as 
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dL = y pac 




/U 


7^+0 +0 

Ut P . . 


dr 


= 2 


- yfiR3 cos ip j dr 


<0. 


-R 


+ y pac 


U?h. 
t 1 


0 +0 sin ijj + 0 cos li; + — 0, 

o s ^ c ^ R 1 


+ ( 


This equation can now be simplified to 


, h. dL = F. (if/> + Y. m., W + X. Q. . (ij;) + X. P. . (ip) 
' 1 1 h ih j j ij 


The generalized force is 


'R 


C'py = y pac 


Ut,hy. 


0 


AfiR ■+ ”00^ '+ 0g sin \p :+ 0^. coS' ip- + 


The generalized coupled torsional stiffness is 

“ih = T f "t \ ®h <'“■ 

The generalized bending stiffness is 

'R 


Qij W = - 2 P^^^ 


h^ y®3^. cos [p dr 


0 


The generalized bending damping is 


'R 


Piy Ci^) = J pac- 


U. h. h. dr 

t 1 J- 
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(A16)' 

\ ■ 

dr 

/ 

(A17) 


y0^Ojdr.' 

(A18) 

'CA19) 


(A20) 


(A2T) 


Blade Torsional Motion : - Considering only two kinds o£ aerodynamic 
feathering moments, and Mp 2 » ^ two-dimensional airfoil [4], the tor- 
sional equations of motion were written as follows: 


a) The first feathering moment, Mp^^, will .occur due to the steady 

aerodynamic lift and the position of the aerodynamic center in 
forward and reverse flow regions. The lift generated by the 
reversed flow can be written as 

dLp = - y ca'p uj dr (A22) 

The first feathering moment can then be written as 

dM„, = ^ c^Ka'p a dr (A23) 

FI 2 t r 

Using the principle of virtual work, the generalized moments 
at each mode can be written as 


'R 


«i ‘"fi. 


12 2 — 
j c^Ka'p 0^ 


U 




dr 


iAa'p J /e. 

+ a. uj 

1 t 


X 'or + cos Ip 


(A24) 


0 : + 0 
P 


dr 


Eq. (A24) can be simplified as 


0i dMpp = G. (Ip) + (iP) + Xj^ (ip) + Y. q. . (ip) 


(A25) 


b) For the second feathering moment, Mp 2 > only one unsteady aero- ' 

dynamic effect is considered and that is the moment due to forces 
caused by the damping of the pitching motion. Assuming that the 
feathering motion occurs about the 1/4 chord-width, the lift due 
to blade pitching velocity has two components: one at 1/4 chord- 

width and one at 3/4 chord-width. The lift force can be expressed 
as 
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(A26) 


dL 


1 2 
7 T c 


ap 


a 

r 


dr 


where 


a = 0+0 
r p 

The lift force in Eq. (A26) generates a feathering moment 


‘^F2 = - T6 ^t “r 

which is valid for the forward flow. Using the principle of 
virtual work, the generalized moments for this case can be 
written as 


'R 


0. dM„_ 
1 F2 


^ c^ ap U. 0 . a dr 
16 ^ t 1 r 


'R 


13^ 

-I6" 


Y. U^ 0. 0. + 0. 0 

3 t 1 j tip 


dr 


Y. p. . + Cn. e 
J ^11 01 P 


Combining Eqs. (A25) and (A29), the generalized torsional moment 
tion can be written as 


6i dMp = Y. p,.(« * Cg,(« * G.w . X, R,,(« * S^. 


») 


+ Y. q. . (4^) 
1 ^11 

where the generalized torsional damping is 

-R 


P • • C4') = - I ap U. 0 . 0 . dr 

^13 I . . 16 ^ t 1 3 

0 


The generalized torsional stiffness is 


(A27) 


(A28) 


(A29) 


equa- 


(A30) 


(A31) 
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qy C« = 


1 c^Ka' p 0 . 0 . dr 

2 ^ t 1 j 


0 


The generalized coupled bending stiffness is 

-R 

12 

c K3. * ^ 


Rj^i(4^) = - I i c^Ka' pyfiR3,, 0^ cos ijj dr 


0 

The generalized coupled bending damping is 

'R 


SkiCW = I l-c^Ka'p hk 0. dr 
0 


The generalized control moment damping is 

-R 


= - I I6 “t *>i 
0 


The generalized control moment is 

'R 


G. W = 


0 


^c^Ka'p U 0. I 0 0 + ASIR 1 dr 

2 t 1 \ t p j 


(A3 2) 


(A3 3) 


(A34) 


(A35) 


(A36) 
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TABLE la. Generalized Structural Mass .0254 m. for Uncoupled Modes 

of CH-34 Blade 
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TABLE lb. Generalized Structural Stiffness (Y^ . = K £.), X 0 = .0254 m. for Uncoupled 
Modes of CH-34 Blade 
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TABLE 2. Uncoupled Generalized Structural Mass , Critical Damping (3.), and Stiffness (y.) 
of CH-34 Blade 
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Fig. 2 Variation of Aerodynamic Bending Stiffness with Blade Azimuth 
Positions, (y = 0.22, k = 0.5, Xq = 0.0). 
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Fig. 3 Variation of Aerodynamic Bending Damping with Blade Azimuth 
Positions, (y = 0.22, K = 0.5, x„ = 0.0). 
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Fig. 4 Variation of Aerodynamic Bending Coupling with Blade Azimuth 
Positions, (p = 0.22, k = 0.5, Xq = 0.0) . 
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Fig. 5 Aerodynamic Stiffness in First Bending Mode at Various Advance 
Ratios, (K = 0.5, Xq = 0.0). 
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Fig. 6 Aerodynamic Damping in First Bending Mode at Various Advance Ratios 
(K = 0.5, Xq = 0.0) . 
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Fig.; 7' Aerodynamic Stiffness Coupling, of First Torsional and First Bending 
Modes at Various Advance Ratios, (k = 0.5, Xg = 0.0). 
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Fig. 8 Aerodynamic Stiffness ,for First Torsional Mode with Various Advance 
Ratios, (K = 0.5, Xg = 0.0). 
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Fig. 9 Aerodynamic Damping for First Torsional Mode with Various Advance 
Ratios, (K = 0.5, Xg = 0.0). 
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Fig. 10 Aerodynamic Stiffness of First Bending and First Torsional Modes at 
Various Advance Ratios, (k = 0. 5, Xg = 0. 0) . 
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Fig. 11 Aerodynamic Damping Coupling' of First Bending and First Torsional 
Modes at Various Advance Ratios, (k = 0.5, Xq = 0.0). 
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Fig. 16 Variation of Coupling Ratios for the Bending Mode with Blade 
Azimuth Positions for Various Eccentricities between Elastic 
and Mass Axes, (u = 0.22, K = 0.5) . 
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Fig. 17 Variation of Coupling Ratios for the Torsional Mode with Blade 
•ivAzimuth Positions for Various Eccentricities between Elastic 
;and Mass Axes, .. (p = 0.22, k = 0.5). 
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18 Variation of Natural Bending Frequency with Blade Azimuth 

Positions for Various Eccentricities between Elastic and Aero 
dynamic Axes, (u = 0.22, Xg = 0.0). 
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Fig. 19 Variation of Coupling Ratios for the Bending Mode with Blade 
Azimuth Positions for Various Eccentricities between Elastic 
and Aerodynamic Axes, (p = 0.22, Xg = 0). 
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Variation of Natural Torsional Frequency with Blade Azimuth 
Positions for Various Eccentricities between Elastic and Aero 
dynamic Axes, (y = 0. 22, Xg = 0. 0) . 
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Fig. 21 Variation of Coupling Ratios for the Torsional Mode with Blade 
Azimuth Positions for Various Eccentricities between Elastic 
and Aerodynamic Axes, (y = 0.22, Xg = 0.0). 
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Fig. 22 Variation of Natural Bending Frequency with Blade Azimuth 

Positions for Various Tip-Speed Ratios, (K = 0.5, Xg =. 0.0). 
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Fig. 27 Variation o£ Coupling Ratios for the Bending Mode with Blade 

Azimuth Positions for Various Tip-Speed Ratios, (k = 0.6, Xg = 0) . 
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Fig. 28 Variation of Natural Torsional Frequency with Blade Azimuth 
Positions- for Various Tip-Speed Ratios, (k =0.5, Xq = 0.0). 
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Fig. 29 Variation of Coupling Ratios for Torsional Mode with Blade 

Azimuth Positions for Various Tip-Speed Ratios, (k = 0.5, Xg = 0.0). 
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Fig. 33 Variation of Coupling Ratios for the Torsional Mode with Blade 
Azimuth Positions for Various Tip-Speed Ratios, (k = 0.6, x-, = 
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Fig. 34 ■ Variation of Natural Bending and Torsional Frequencies with Blade 
Azimuth Positions for Various Eccentricities between Elastic and 
Mass Axes, (y = 0.22, k = 0.4). 
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Fig, 35 Variation of Coupling Ratios for the Bending Mode with Blade 
Azimuth Positions for Various Eccentricities between Elastic 
and Mass Axes, (p = 0.22, K = 0.4). 
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